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Abstract 

In environmental sciences, it is often of interest to assess whether the dependence between 
extreme measurements has changed during the observation period. The aim of this work 
is to propose a statistical test that is particularly sensitive to such changes. The resulting 
procedure is also extended to allow the detection of changes in the extreme-value dependence 
under the presence of known breaks in the marginal distributions. Simulations are carried 
out to study the finite-sample behavior of both versions of the proposed test. Illustrations 
on hydrological data sets conclude the work. 
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1 Introduction 

The study of extremes is of importance in many environmental applications. Prominent exam¬ 
ples are the analysis of floods (Hosking and Wallis, 2005), heavy rainfalls (Cooley et ah, 2007) 
and extreme temperatures (Katz and Brown, 1992). Many of these problems are intrinsically 
multivariate; for instance, the severity of a flood depends not only on its peak flow, which is 
considered in many univariate flood studies, but also on its volume and its duration (Yue et ah, 
1999). Catastrophic flood events typically occur when more than one of these variables is taking 
a high value and therefore, the analysis of the joint behavior is of key importance. In a river 
system, where flood data are collected from a number of stations, inference at a specific loca¬ 
tion can be greatly improved by incorporating observations from neighboring stations (Hosking 
and Wallis, 2005). Similarly, extreme temperatures are commonly studied at several stations 
simultaneously. 

In most of these environmental applications, it is common practice to assess the extreme 
observations by (modifications of) the annual block maxima method, popularized in the classi¬ 
cal monograph by Gumbel (1958). Univariate observations collected on, say, a daily basis are 
aggregated by taking maxima over a longer time period, usually a year for stationarity consider¬ 
ations, resulting in a sample of maxima capturing most of the extreme outcomes. Thanks to the 
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extremal types theorem, univariate block maxima are approximately distributed according to 
the (parametric) generalized extreme distribution, while the dependence within a vector of block 
maxima can be described approximately by a (nonparametric) extreme-value copula (de Haan 
and Ferreira, 2006). The resulting joint distributions are called multivariate extreme-value 
distributions and serve as a widely accepted model for multivariate block maxima. 

In statistical applications, it is common practice to assume that the time series of block 
maxima is stationary, temporarily independent and that its stationary distribution is exactly 
of the extreme-value type. Except for stationarity, the above assumptions are theoretically 
justihed by the results of Biicher and Segers (2014) if the focus is solely on the dependence 
structure. In the current paper, we address the issue of stationarity by developing tests for 
change-point detection within the multivariate contemporary distribution. 

More precisely, assuming that we observe a sample of independent multivariate observations 
Xi,...,X„, where each Xi follows a multivariate extreme-value distribution whose c.d.f. is 
denoted , we develop a test for the hypothesis 

Ho : = ... = (1.1) 

against alternatives involving the non-constancy of the c.d.f.s. Since the univariate version of 
this problem has been treated, for instance, in Jaruskova and Rencova (2008) using results from 
Chapter 1 of Csorgd and Horvath (1997), we will be particularly interested in the multivariate 
setting throughout this paper. 

Outside of the extreme-value framework, there is a huge amount of literature on detecting 
deviations from Ho, see the classical monograph of Csorgd and Horvath (1997) or Aue and 
Horvath (2013) for a recent review. It is useful to note that, by Sklar’s theorem (Sklar, 1959), 
we can rewrite Ho as Ho,m 77o,c, where 

Ho,m : ..., have time-homogeneous marginal c.d.f.s, (1-2) 

Ho,c : ..., have the same copula (i.e., dependence). (1-3) 

Roughly speaking, tests for Ho can be divided into two groups: tests that are powerful against 
deviations from Ho,m but which are rather insensitive to deviations from Ho,c, and vice versa, see 
Biicher et al. (2014) for a discussion. In the present setting of observing data from a multivariate 
extreme-value distribution, the tests considered for instance in Jaruskova and Rencova (2008) 
can be used to detect deviations from Ho,m, whence it will be our main concern to design a 
test that is particularly sensitive to deviations from Ho,c when C is known to be an extreme- 
value copula. Note that none of the existing tests for changes in the copula (see, e.g., Bucher 
et ah, 2014; Kojadinovic et ah, 2015, among others) incorporates the latter information, whence 
an improvement in the power properties seems possible. In fact, our simulation study partially 
reported in Section 5 suggests that our proposed testing procedure is indeed superior to existing 
methods. 

The test tailored to deal with extreme-value dependence that we propose is however affected 
by the same limitations as the aforementioned more general testing procedures: it can be used 
to reject Ho,c only if Ho,m holds. In some situations, although there are reasons to consider 
that Ho,m is not true, it is still of interest to assess whether Ho,c holds or not. For instance, in 
the hydrological applications to be presented in Section 6, the construction of dams during the 
observation period suggests that there might be potential breaks in the marginal distributions 
of extreme peak flows or volumes, while it is still of interest to assess whether the dependence 
between the variables of interest has changed or not. A second contribution of this work is thus 
to propose an extension of the studied testing procedure that can be used to detect deviations 
from Ho,c under certain types of simple departures from Ho,m- 
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The remaining part of the paper is organized as follows: The second section is devoted 
to mathematical preliminaries about extreme-value copulas. In Section 3, a first version of 
the testing procedure is presented along with theoretical results establishing its asymptotic 
validity under Hq. An extension of the studied test that can be used to detect change in the 
extreme-value dependence under known marginal breaks is proposed in Section 4 along with 
generalizations of the theoretical results of Section 3. Section 5 partially reports the results 
of extensive Monte Carlo experiments. Two illustrations on hydrological data sets are finally 
presented in Section 6. All proofs are deferred to a sequence of appendices. 


2 Notation and mathematical preliminaries 

In this section, we set the notation and gather some mathematical preliminaries needed through¬ 
out the paper. Let X = {X\,... ^ Xd) be a d-dimensional random vector with continuous 
marginal c.d.f.s Ti,...,Trf. By Sklar’s theorem (Sklar, 1959), the joint c.d.f. H oi X can be 
uniquely decomposed as 

H{x) = C{Fi{xi),...,Fd{xd)}, a; = (xi,... ,Xrf) 6 (2.1) 

where the so-called copula C : [0,1]'^ ^ [0; 1] is the joint c.d.f. oi U = {Ui,... ,Ud) with 
Uj = F,{Xj). 

Throughout the paper, we will assume that the copula C in (2.1) is an extreme-value copula. 
As a consequence of the results in Pickands (1981), these copulas can be parametrized by a 
function A : Sd-i [1/d, 1] on the (d — l)-dimensional unit simplex Sd-i = {t = {t 2 , ..., t^) e 
[0,1]'^“^: t 2 + ■.. + td ^ 1}, usually referred to as the Pickands dependence function. More 
precisely, we have 


C{u) = exp 



A 


log U 2 

SjlllogUj 


log Ud 

’ 2^=1 log 


( 2 . 2 ) 


for any u e (0,1]‘^\{(1,..., 1)}. If relation (2.2) is met, then A is necessarily convex and satisfies 
the boundary condition max{l — Xij= 2 ^ 2 , • • •,trf} ^ A{t) ^ 1 for all t = {t 2 ,... ,td) ^ •Sd-i- 
The latter two conditions are, however, not sufficient to characterize the class of Pickands 
dependence functions unless d = 2, see, e.g., Beirlant et al. (2004) for a counterexample and 
Ressel (2013) for recent results concerning the case d > 2. 

Assuming to observe a sample Xi, i = 1,... ,n, of serially independent random vectors such 
that Xi has copula and corresponding Pickands dependence function A^^\ we can use the 
representation in (2.2) to rewrite Pfo,c in (1-3) equivalently as 

Ho,c : = ... = A^^^ = A. (2.3) 


The test statistics in the subsequent sections will be particularly designed for detecting devia¬ 
tions from this hypothesis. 

In the rest of the paper, for notational convenience, we will work mostly in dimension d = 2. 
Furthermore, given a set T, (£°°(T), || • loo) will denote the space of real-valued, bounded functions 
on T equipped with the supremum norm || • ||oo. The arrow denotes weak convergence in 
the sense of Hoffmann-Jprgensen, see van der Vaart and Wellner (1996). 
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3 Test statistics for d = 2 under stationarity of the marginals 


We begin by restricting ourselves to the case of dimension d = 2 and assume to observe an 
independent sample (W, y), z = 1,..., n, such that (W, y) has unknown c.d.f. copula 
and continuous margins and , respectively, where each is assumed to be an extreme- 
value copula with Pickands dependence function : [0,1] ^ [1/2,1]. 

We aim at developing tests for change-point detection that are consistent against deviations 
form in (2.3) provided ?^o,m in (1-2) holds. 

Our test statistic will be of the CUSUM-type and will be based on the rank-based estimator 
of A proposed by Ferreira (2013). The underlying idea of that estimator is as follows; if ([/, V) 
is distributed according to some extreme-value copula C with Pickands dependence function A, 
then 

E{max(t/lV--).rV‘)} 

' 1-lE{max((/l/(l-<),yl/<)}’ 

with the convention that = 0 for any u e (0,1). The function 


S{t) = E I max |, t 6 [0,1], 


(3.1) 


is simply an expected value, whence defining an estimator for S under the null hypothesis 
Tio = T-Lo^rn ^o,c in (1-1) is straightforward once we have estimated the unobservable sample 
(C/j, Vi) = (F(W)) G{Yi)), i = 1,... ,n. To do so, it is common to compute the scaled ranks 


U -| it 

Ul:n,i = —— 1 ^ ^ 1 (Tj ^ Tj) : i = 

” + S = 1 ” + i = l 


(3.2) 


frequently referred to as pseudo-observations from the unknown copula C. Then, a natural 
estimator of S is simply 

1 ^ 

Sl:nit) = - 2 ^ ^ [ 0 , !]• 

2=1 

The corresponding estimator of A, namely 

il:n(t) = Si..nit)/{1 - 5i;„(t)}, t 6 [0, 1], 


was shown to be consistent and was investigated empirically in Ferreira (2013). As a by-product 
of our work, we establish the asymptotic distribution of the process ^/n{Al■n—A) in the following 
proposition proved in Appendix A. To the best of our knowledge, this result is new and might 
be of independent interest, for instance for the construction of uniform confidence bands. 

Proposition 3.1. Suppose that (Aj,y), i = l,...,n, are i.i.d. from a bivariate distribution 
with continuous marginal c.d.f.s and extreme-value copula C whose Pickands dependence func¬ 
tion A is continuously differentiable on (0,1). Then, in the normed space (£°°([0,1]), jj • |loo), 
^/n{Al■n - A) -> Lc, where 

^c{t) = -{1 + Mdf f Cc(0,1, yd dy 

Jo 

and Cc'(0,1,-,-) is the weak limit of the empirical copula process (see, e.g., Segers, 2012) as 
defined in Proposition 3.2 below. 
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In order to define a CUSUM-type procedure for testing Hq in (1.1), we consider the following 
subsample analogues of Si-n based on subsequences Y),),..., (X^, Y)), 1 ^ A; ^ ^ n, 

namely 

= i-k + l 2 (3.3) 

i=k 

where 

i £ 

j=fc j=k 

with the convention that Sk-i = 0 for k > 1. The corresponding subsample estimators of A are 
then simply 

Aki{t) = Sk-,i{t)/{l - Sk-.t{t)}, t e [0,1]. (3.5) 

Under T-Lq, the difference between Ai-k and Ak+i-n should be small for any A; = 1,..., n — 1, 
which suggests to base test statistics for "Hq on the process 

^n{s,t) = {A:M(A) “ \ns\ + l:n{t)} , (3.6) 

for {s,t) e [0,1]^. Typical test statistics would be given by Kolmogorov-Smirnov or L^-type 
functionals of D^. Throughout this paper, we focus on the hybrid version 

Sn,A = max I {Bn(k/n,t)}^ d/i(t), (3.7) 

l^k<n Jj-Q 1 ] 

where /r denotes some finite measure on [0,1]. In the finite-sample experiments of Section 5, 
we use /i = T~^ Xiter some finite grid T = {ii,..., ir} c [0,1], where 6t is the Dirac mass 
at t. A related two-sample statistic for detecting breaks at some pre-specified time point k* is 
simply given by 


5n,A(fe*)= {B^ik*/n,t)f (3.8) 

J[o,i] 

Note that the aforementioned test statistics do not incorporate the information that the marginal 
distributions are (or should be close to) generalized extreme-value (GEV) distributions. There 
are several reasons for ignoring that information: First of all, since we are only interested in 
the dependence, it seems natural to ignore any marginal information. This is theoretically 
justified by results in Genest and Segers (2010), where it is shown that rank-based estimation 
of copulas may be substantially more efficient than estimation based on even exact knowledge 
of the marginal distributions. Second, even if we opt for a semiparametric estimation of the 
dependence structure based on a parametric estimation of the marginal distributions, we could 
expect the estimation of the copula to be quite negatively affected by the fact that estimators 
of GEV parameters have frequently a large variance. Last but not least, extreme-value copulas 
are also of interest outside of the genuine extreme-value framework, where margins are not 
necessarily of the GEV-type. 

The following proposition, proved in Appendix A, establishes weak convergence of the key 
process On in (3.6) under T-Lq in (1.1). It is essential for deriving the weak limit of the preceding 
test statistics. 
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Proposition 3.2. Suppose that %q holds and that A is continuously differentiable on (0,1). 
Then, in the normed space 1]^), || • |loo), On Oc, where 

Oc(s, t) = {1 + A{t)f r sCc{s, l,y^~\y^) - (1 - s)Cc(0, s, y^~\ y*) dy. (3.9) 
Jo 

Here, Cc denotes a centered Gaussian process on Ax [0, 1]^, with A = {(s, s') e [0,1]^ : s ^ s'}, 
defined through 

Cc{s, s', u, v) = {IBc('s', u, v) - IBc('S, u,v)} - Cffu, 'u){]Bc(s', u, 1) - Bc(s, u, 1)} 

- C 2 {u,v){Mc{s', l,v) -]Bc('S, 1,^^)} 

if {u,v) 6 (0,1)^ and Cc{s,s',u,v) = 0 else, where Be denotes a tight, centered Gaussian 
process on (£“([0,1]^), || • ||oo) with covariance kernel 

E{Be(s, u, v)Mc{s', u , u')} = min(s, s')[C{min(u, n'), min(t!, t?')} — C{u, v)C{u', i;')], 

and where Cj, j = 1, 2, denotes the jth first-order partial derivative of C. 


The limit process Oe in (3.9) depends in an intractable way on the unknown copula C, and, 
as a consequence, so do the limit distributions of the test statistics Sn,A in (3.7) and Sn,A{h*) 
in (3.8). The following multiplier bootstrap approximations, initially proposed by Scaillet (2005) 
and Remillard and Scaillet (2009) in a copula setting, allow for the derivation of suitable ap¬ 
proximations of the critical values. 

Let B be some large integer and let i = 1,... ,n, b = 1,... ,B, denote i.i.d. standard 
normal random variables. Motivated by recent results in Bucher and Kojadinovic (2014) and 
Bucher et al. (2014) on the multiplier bootstrap in a sequential setting, we can approximate the 
Gaussian process Be appearing in the limit De defined in (3.9) by 


[ns'J 


,S,U,v) ypn ^ ^ f^[nsJ + l:[ns'J (^) ^) !■) 




(3.10) 


for iys,s',u,v) 6 A X [0,1]^, where, for 1 ^ A: ^ ^ n, C'fc:£ denotes the empirical copula based 

on the sample [X^, Yff), ■ ■ ■, {X^, Yfij, that is. 


Ck-.l{u, = {Ukl,i ^ u, Vk:i,i ^ u), (u,u) 6 [0, 1]^, (3.11) 

i=k 

with the convention that Ck-i = 0 ii k > £. If we replace Be by B^f^ and all other unknown 
quantities in the definition of De by natural estimators, then some standard calculations, carried 
out explicitly in the proof of the next proposition, suggest to define bootstrap replicates of 
as 



{1 + Ai-n{t)}'^ X 




(3.12) 
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for (s, t) e [0,1]^, where, for any 1 ^ A: ^ ^ ^ n, 


Wk:e,i{t) = mk-i{t) - mk-.e,i{t) + {uk-.e,i{t) - Uk-,e{t)} ^^'^^^l + {vki,i{t) - (3.13) 

bk:i{t) dk-i{t) 




with mk-j,Uk-.£ and Vk-.i denoting the arithmetic mean (over i = k, oi 

rhkMt) = 

respectively, and with 

dki{t) = Ak-iit) - bk,i{t) = Ak,i{t) + t, 

Ck-£{t) = Ak-At) + (1 - t)i^„(t), dk,i{t) = Ak,i{t) + 1 - t, 

where, for some positive sequence /in i 0 such that infn^i hny/n > 0, 

A’k-.e,nit) = min[max{^'fc,^ ,,(t), -1}, 1], te [0,1], (3.14) 

with 

Ak-.e,ni^) ~ |^A::£(/ + hn) — Ak:£{t — /ln)| , (3.15) 

for t 6 {hn, 1 - hn), while A'j..^{b) = A'f..^ A^n) for t ^ hn and A'^.g^ nit) = “ ^n) for 

t ^ 1 — hn- 

The following proposition, proved in Appendix B, establishes the asymptotic validity of the 
above resampling scheme under T-Lq in (1.1). 

Proposition 3.3. Under the conditions of Proposition 3.2, 

in (^”([0,1]^), II • lloo)'®’^^, where denote independent copies o/Dc- 


Bootstrap analogues of the test statistics in (3.7) and (3.8) can be defined by replacing 
On by O),^) in the corresponding definitions. Through the continuous mapping theorem, we 
immediately obtain that, under T/q, the random vector {Sn,A, ..., e weakly 

converges to a vector with i.i.d. components, each component having the same distribution 
as sup^g|-Q y 5[o 1 ] ^)'^/^(^)- Hence, a test rejecting 77o at the significance level a if Sn,A 

exceeds the [(1 — a)B\-ih. order statistic of S!^\,..., asymptotically keeps its level for 
n ^ 00 followed by i? ^ oo (see Appendix F in Biicher and Kojadinovic, 2014). 

At the cost of a more complex notation, the results presented in this section can be extended 
to the case d > 2. The main steps are given in Appendix C. 


4 Test statistics for d = 2 under known marginal breaks 

It is only if ?7o,m in (1-2) holds that the tests developed in the previous section can be used to 
reject ?7o,c in (2.3). Empirical evidence of the latter fact will be given in Section 5. In some 
applications, such as those presented in Section 6, there might be reasons to believe that (1.2) 
does not hold. The aim of this section is to propose an adaptation of the tests derived in 
the previous section that will be consistent against ?7o,c in (2.3) when there are known abrupt 
changes in the margins. 
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Recall that we assume to observe an independent sequence {Xi,Yi), i = where 

{Xi,Yi) has copula and continuous marginal c.d.f.s and respectively. Before 

considering a more general framework, we shall first focus, for pedagogical reasons, on the 
following simple hypothesis for the margins: 


n 


l,m 


there exists a known 0 e (0,1) such that 
i?(i) = ... = i?(h^'l) ^ _p(h^'l+i) = ... = F^'^\ 
qO-) = ... = (^(R^'l) ^ G'lR^'l+i) = ... = G^'^\ 


(4.1) 


that is, both margins change abruptly at the same time point \n9\. As we continue, we shall 
also use the notation m = \n0\. In the illustrative applications presented in Section 6, this time 
point corresponds to the construction of a dam on a river. 

In order to adapt the testing procedures derived in the previous section to the above setting, 
we need to propose an appropriate version of the test statistics Sn,A in (3.7) and Sn,A{k*) in (3.8) 
under 7fi,mn?7o,cj where ?7o,c is defined in (2.3). Mimicking the developments carried out in the 
previous section, it is natural to start from adequate versions of the estimators in (3.3). Hence, 
assume that l-Li^rn ’Ro,c holds, and let (Aj, Yi), i = k, ...,£, he some subsample. If k < m < i, 
a natural estimator for S in (3.1) is given by a convex combination of the estimators defined 
in (3.3), one based on the subsample from k to m, the other on the subsample from m + I to £. 
In other words, meaningful analogues of the estimators in (3.3) are, for any t e [0,1], 


_ j + |z:^*gm+l:^(t), if fc < m < £, 


siAt) = 


SkAt), 


otherwise. 


(4.2) 


Proceeding as in the previous section, the corresponding subsample estimators of A are then 
simply 




sUt) 


te [0,1]. 


(4.3) 


1 - siAt) 

It is easy to verify that the formulas in (4.2) and (4.3) coincide with those in (3.3) and (3.5), 
respectively, provided that the pseudo-observations in (3.4) are replaced by appropriates ones 
taking into account the break in the margins at time point m, namely, for i = k,... 


Ue,k.i,i = " 


Uk:£,i 

Uk-.m,i 

Um+l:i,i 


\im^ {fc, . . . £}, 

if m e {fc, ... A} s-iid i ^ m, 

if m e {k,... A} i > m, 


(4.4) 


and similarly for V 0 ^k:t,i- Thus, from a practical perspective, once the above adapted pseudo¬ 
observations are computed, the computer code for the simpler setting considered in the previous 
section can be fully reused. 

It follows that natural generalizations of the process and the statistics Sn,A and Sn,A{k*), 
defined in (3.6), (3.7) and (3.8), respectively, are simply 

<(0,i) - . (».*) <5 [0,1]^ 

and 

= max f {D^(/c/n,t)j d/i(t), SiAk*)={ |D®(fe7n,t)j d/x(t), (4.5) 

respectively. 

The following proposition, proved in Appendix A, establishes the limit distribution of 
under T-Lo,c where 'Ho,c and are defined in (2.3) and (4.1), respectively. 
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Proposition 4.1. Assume that ?^o,c holds and that the Pickands dependence func¬ 

tion A associated with C is continuously differentiable on (0,1). Then, in the normed space 
1]^), II • j|oo)j where Dc is defined in (3.9). 

Hence, we see that the limit distribution of the process D® under TLo,c ^ 'Hi,m does not 
depend on the marginal break point 9 and coincides with that of the process in (3.6) whose 
asymptotics were studied in Proposition 3.2 under Ho,c TLo,m, where TLo^m is defined in (1.2). 

To carry out the tests based on S'® ^ and we need bootstrap replicates of the pro¬ 
cess D®. The latter are denoted by b = 1,... ,B, and are defined exactly as in (3.12), 

except for the underlying pseudo-observations which are computed as in (4.4). The next result, 
proved in Appendix B, establishes the asymptotic validity of the proposed resampling scheme 
under Pfo,c Tdi,m- 

Proposition 4.2. Under the conditions of Proposition f.l, in (£“([0,1]^), || • ||oo)^''’^, 

where Dc is defined in (3.9) and denote independent copies o/Dc. 


Multiplier bootstrap replicates of the test statistics in (4.5) are defined mutatis mutandis, the 
only difference with the formulas of Section 3 being again the pseudo-observations which are 
computed using (4.4). As previously, the null hypothesis is rejected at the significance level a 
if 5®^ (resp. S^^{k*)) exceeds the [(1 — Q;)Hj-th order statistic computed from the sample of 
its B multiplier bootstrap replicates. 

We end this section by three remarks on further possible extensions. 


Remark 4.3. (More than one marginal break) In the preceding parts of this section, 
we restricted ourselves to the case of exactly one abrupt break that affects both marginal 
distributions. At the cost of a more complex notation, the previous results can all be simply 
extended to a more general setting. For i? e N, i? ^ 1, let 

©i? = {(6'o, 9i,. .. ,6r+i) 6 : 0 = 6q < 6i < ■ ■ ■ < Or < 6r+i = 1}. 


A vector 6 e Qr should be interpreted as encoding (rescaled) time points where at least one of 
the marginal distributions changes abruptly. Specifically, let 


r 


n 


(R) . 

l,m * 


< 


there exists a known 0 e Qr such that, for any r = 0,R, 
p{[ner\ + l) = _ = p{[n9r+i\) 

Q{[ner\ + l) ^ ^ Q{[ner + l\) 

and such that, for at least one of the margins, 
the c.d.f.s at and + 1 are different. 


Identifying the parameter 6 e (0,1) used previously in this section with the vector 6 = (0, 9, 1) e 
01, we see that in (4.1) implies above. To derive an extension of the test based 

on 5® ^ in (4.5) that can be used to reject Ffo,c in (2.3) under above, it merely suffices 
to adapt the definition of the pseudo-observations in (4.4) to this more complex situation. For 
that purpose, let = [n9r\ for all r = 0,..., R -I- 1, and, for any f = 1,..., n, let 

tUa = m-A (i) = max{mr. e {mo ,..., m/j} : < i}, 

ruv = ruv (i) = min{mr. 6 {mi,..., m/j+i} : mr ^ i]. 
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such that i e {m* + 1,... ,mv}. Then, analogously to (4.4), we can define pseudo-observations 
adapted to for 1 ^ A: ^ ^ n and i = k,... ,i, as 


Ue,ki,i 


Uk:£,i 

Uk-.rtiy ,i 
^rriA +l:mv ,i 


if < k and i < , 

if ^ k and (. < , 

if < k and i ^ , 

if ^ k and £ ^ ruv, 


and similarly for Ve^k-.i,i- All the formulas from Section 3 remain identical up to the use of the 
above definition for the pseudo-observations. In addition, the proofs of the theoretical results 
stated in this section extend easily but are more cumbersome to write. For the sake of brevity, 
we omit further details. 


Remark 4.4. (Extension to PQD copulas) The previously studied tests could be extended 
to copulas that are not necessary of the extreme-value type, for instance to positive quad¬ 
rant dependent (PQD) copulas (see, e.g., Nelsen, 2006, Chapter 5), that is, copulas satisfying 
C{u,v) ^ uv for all {u,v) 6 [0,1]^. Starting from (2.2), it can be easily verified that extreme- 
value copulas are PQD. The extension of the tests to PQD copulas relies on the fact (used in 
the proofs given in Appendices A and B) that the Pickands dependence function associated 
with an extreme-value copula C can be expressed as 


A{t) 


SlC{y^-\y^)dy 


te [0,1]. 


The previous functional of C remains well-defined outside of the extreme-value framework as 
long as 5o C{y^~^,y^)dy A 0. The PQD condition, for instance, ensures that C{y^~^,y^) ^ y for 
all {y,t) £ (0) 1) ^ [0) l]j and therefore the existence of A. 

The only change necessary to extend the previously studied tests to PQD copulas concerns 
the estimation of the partial derivatives Ci and C2 of the copula C that is required for car¬ 
rying out the multiplier resampling scheme. The approach adopted in Section 3 (see (3.15)) 
could be replaced by the more classical one consisting of using finite-differences based on the 
empirical copula. This would however lead to significantly less convenient formulas as far as 
implementation is concerned. Given that the meaning of the functional A is unclear outside of 
the extreme-value framework and that tests for change-point detection for non extreme-value 
copulas already exist (see, e.g., Bucher et ah, 2014; Kojadinovic et ah, 2015), we do not pursue 
this further. 


Remark 4.5. (Estimation of the marginal change-points) An even more general frame¬ 
work regarding marginal change-points is to assume that their number is known but not their 
position. For instance, under the assumption of one unknown marginal (scaled) change-point 
9 e (0,1) in the first margin, its value could be estimated under some additional assumption on 
the nature of the change-point (change in mean, in variance, etc). A sensible estimator 9n of 6 is 
then given by k/n, where k is the value offce {l,...,n — 1} that maximizes a corresponding nat¬ 
ural max-type change-point statistic. Under the additional assumption that 9n — 9 = Op(l/n) 
(see, e.g., Diimbgen, 1991), we conjecture that the limiting behavior of the statistic is the 
same as that of 5® ^ in (4.5). The theoretical justification seems to be quite involved, and, since 
we do not need this extension for the hydrological applications presented in Section 6, we do 
not pursue this further. 
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5 Simulation study 


Simulations were carried out in order to evaluate the finite-sample performance of the tests 
studied in Sections 3 and 4. For the sake of simplicity, we only focused on the tests based on 
Sn,A in (3.7) and S'® ^ in (4.5), as the finite-sample behavior of the corresponding two-sample 
tests should be strongly related. Recall that the tests based on Sn,A and S® ^ are procedures for 
testing Hq in (1.1) designed to be particularly sensitive to departures from ?7o,c in (2.3). The 
former (resp. latter) should not be used to reject ?7o,c unless 77o,m in (1.2) (resp. in (4.1)) 
holds. 

The rejection rates of the two tests were estimated from samples drawn from bivariate dis¬ 
tributions whose copulas are of the form 

Ca,^{u,v) = {u,v) e [0,1]^, (4.6) 

where is a symmetric extreme-value copula with parameter e M and a = (oi, 02 ) e [0,1]^ 
is a parameter controlling the amount of asymmetry. The above copula construction principle 
is frequently referred to as Khoudraji’s device (Khoudraji, 1995; Genest et ah, 1998; Liebscher, 
2008). As long as is an extreme-value copula, so is its potentially asymmetric version 

Given that there is hardly any practical difference among the existing bivariate symmetric 
parametric families of extreme-value copulas (see Genest et ah, 2011, for more details), 
in (4.6) was taken to be the Gumbel-Hougaard copula with parameter d. 

The finite-sample performance of the tests based on Sn,A and S'® ^ was compared with that 
of two other tests for T-Lq designed to be particularly sensitive to ?7o,c in (1.3): the test based 
on the empirical copula studied in Bucher et ah (2014) (statistic Sn) and the test based on 
Spearman’s rho considered in Kojadinovic et ah (2015) (statistic Sn,i)- Both tests rely on 
multiplier bootstraps for the computation of approximate p-values. They will be referred to as 
the test based on Sn,c and the test based on Sn,p, respectively. These procedures however do 
not assume the underlying dependence structures to be of the extreme-value type. The former 
is sensitive to all kind of changes in the underlying copula, while the latter is only sensitive to 
changes in Spearman’s rho. As for the test based on Sn,A, they should not be used to reject 
?7o,c unless ?7o,m in (1.2) holds. 

All the tests considered in our numerical experiments were carried out at the 5% signihcance 
level using B = 1000 multiplier bootstrap replicates. The values 50, 100 and 200 were considered 
for the sample size n. The measure fj, involved in the definition of Sn,A in (3.7) and S'® ^ in (4.5) 
was taken equal to 9“^ hj/iQ. The bandwidth hn in (3.15) was set to lO^^/yhr. With 
the illustrations of Section 6 in mind, the values 0.25 and 0.5 were considered for 9. The 
computations were carried out using the R statistical system (R Gore Team, 2015), and the R 
packages copula (Hofert et ah, 2015) and npcp (Kojadinovic, 2014). 

Empirical levels of the tests based on Sn,Aj Sn,c and Sn,p Columns 5-7 of Table 1 report 
the rejection rates of the three tests estimated from 4000 random samples generated under Tio 
from c.d.f. Ca,& in (4.6) for various values of a and d. As one can see, the empirical levels are 
overall reasonably close to the 5% nominal level in all settings for which Kendall’s tau r of Ca,,? 
is strictly smaller than 0.6. For r ^ 0.6, the three tests are overall too conservative, although, 
as expected, the empirical levels improve as n increases. 

Empirical power of the tests based on Sn,A^ Sn,c and Sn,p nnder changes in the 
copula only The right plot of Figure 1 displays the rejection rates of the three tests estimated 
from 2000 samples of size n = 100 generated under 77o,m ~^T~Lo^c such that, for each sample, 
the first (resp. last) 50 observations were drawn from c.d.f. (4.6) with (a, d) = (0,0, 2) (resp. 
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n 

a 

d 

r 

Sn,A 

Sn,C 

Sn,p 

c<0.25 

^n,A 

c-O.S 

^n,A 

50 

(0,0) 

1 

0 

4.9 

6.3 

5.6 

7.6 

4.2 



1.25 

.2 

6.7 

6.2 

5.8 

7.9 

7.0 



1.67 

.4 

5.8 

4.4 

3.9 

6.6 

6.1 



2.5 

.6 

4.0 

3.7 

2.4 

5.6 

4.7 



5 

.8 

3.6 

2.4 

0.9 

8.2 

2.7 


(0,.3) 

4 

.56 

4.5 

4.7 

3.1 

5.2 

5.5 

100 

(0,0) 

1 

0 

5.5 

5.1 

5.8 

7.7 

5.4 



1.25 

.2 

6.3 

5.4 

6.2 

7.4 

6.9 



1.67 

.4 

6.2 

4.3 

4.4 

6.2 

6.6 



2.5 

.6 

5.4 

3.0 

2.9 

6.0 

5.5 



5 

.8 

2.0 

2.2 

1.0 

4.0 

2.6 


(0,.3) 

4 

.56 

4.5 

4.2 

3.8 

4.5 

5.0 

200 

(0,0) 

1 

0 

5.0 

4.3 

4.8 

6.2 

5.6 



1.25 

.2 

6.0 

4.8 

5.8 

6.4 

6.4 



1.67 

.4 

5.9 

4.0 

4.9 

6.4 

6.2 



2.5 

.6 

3.6 

2.8 

3.1 

4.4 

4.4 



5 

.8 

2.6 

1.3 

2.0 

3.4 

3.4 


(0,.3) 

4 

.56 

4.8 

4.0 

4.3 

5.2 

5.2 


Table 1: Rejection rates of 1-Lq in % estimated from 4000 random samples generated under T-Lq 
from c.d.f. Ca,,? in (4.6). The column r gives the value of Kendall’s tan of the copula Ca,^. 



Figure 1: Left: Pickands dependence function of the GumbeKHougaard copula with parameter 
2 + (I'd. Right: Rejection rates of the tests based on Sn,A i°), Sn,c (^) and Sn,p (+) versus 
(I'd estimated from 2000 bivariate samples of size n = 100 such that, for each sample, the first 
(resp. last) 50 observations were drawn from a GumbeKHougaard copula with parameter 2 
(resp. 2 + d-d). 
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0.0 0.2 0.4 0.6 0.8 1.0 



t 


da 


Figure 2: Left: Pickands dependence function associated with the copula in (4.6) with 
a = (max(0.4 — da, 0), max(da — 0.4,0)), for da e {0, 0.4,0.8}, and d set to keep Kendall’s tau of 
Ca,i 9 equal to 0.5. Right: Rejection rates of the tests based on Sn,A i°), Sn,c and Sn,p (+) 
versus da estimated from 2000 bivariate samples of size n = 200 such that, for each sample, the 
hrst (resp. last) 100 observations were drawn the above mentioned copula with da = 0 (resp. 
da e {0, 0.1,..., 0.8}). 


(0,0,2 + dd)). As expected, the test based on Sn,A is more powerful than its two competitors 
in this simple setting. 

To investigate the influence of asymmetry on the power of the three tests, as a second ex¬ 
periment, we considered again the copula Ca,§ in (4.6), but this time with parameter a defined 
as (max(0.4 — da, 0), max(do — 0.4, 0)), for da e {0, 0.1,..., 0.8}, and with parameter d set to 
keep Kendall’s tau of Ca,'& equal to 0.5. The corresponding Pickands dependence functions for 
da 6 {0, 0.4, 0.8} are represented in the left plot of Figure 2. The right plot displays the rejection 
rates of the tests based on Sn,A, Sn,c and Sn,p versus da estimated from 2000 samples of size 
n = 200 such that, for each sample, the first (resp. last) 100 observations were drawn from the 
above mentioned copula with da = 0 (resp. da e {0,0.1,..., 0.8}). Although the rejection rates 
are overall relatively low, the test based Sn,A is by far the best. The fact that the test based on 
Sn,p has no power against such alternatives is due to the fact that Spearman’s remains almost 
constant. 

Empirical power of the tests based on Sn,A^ Sn,c and Sn,p under an abrupt change in 
one margin only Table 2 reports rejection rates of T-Lq estimated from 1000 bivariate samples 
of size n generated under ^o,c where ?^o,c and 'Hi^m are defined in (1.3) and (4.1), 

respectively, such that, for each sample, the first \n0\ (resp. last n — \n9\) observations were 
drawn from a c.d.f. whose copula is the Gumbel-Hougaard, whose first margin is GEV with 
parameters ^ = 20, a = 10 and 7 = 0.25 (resp. // = 20 -I- d/r, a = 10 and 7 = 0.25), and whose 
second margin is standard normal (the results are unaffected by the choice of the second margin 
since the test is rank-based). 

All three tests have little power against such alternatives when the shift d// in the location 
parameter of the first margin is relatively small {dfi = 5). This is a desirable property since 
the tests were designed to be sensitive to departures from 'Ho,c- Higher rejection rates were 
obtained for d/r = 15 and when the dependence is moderate or high, in particular if the (scaled) 
change-point in the first margin is non-central (d = 0.25). The latter results illustrate the fact 
that the procedures based on Sn,A, Sn,c and Sn,p are tests for Ho and that one should not use 
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II 

0 


9 

= 0.25 


dfjL 

n 

T 

Sn,A 

Sn,C 


Sn,A 

Sn,C 

Sn,p 

5 

50 

0 

5.5 

6.2 

6.0 

5.5 

6.3 

6.2 



0.25 

6.6 

6.1 

6.2 

5.2 

4.1 

3.7 



0.5 

4.3 

3.3 

2.5 

4.4 

2.8 

2.7 



0.75 

3.9 

2.2 

1.2 

2.7 

2.5 

1.1 


100 

0 

5.1 

5.1 

5.4 

5.0 

5.9 

6.0 



0.25 

5.1 

3.6 

5.1 

6.7 

4.0 

4.9 



0.5 

5.5 

4.1 

4.3 

4.7 

2.9 

3.3 



0.75 

3.3 

1.0 

0.6 

4.5 

2.8 

1.5 


200 

0 

5.6 

4.8 

5.5 

5.6 

5.8 

5.4 



0.25 

4.8 

5.1 

5.1 

3.9 

3.6 

3.4 



0.5 

4.6 

2.5 

2.9 

5.2 

3.4 

4.0 



0.75 

4.0 

0.8 

1.0 

5.0 

3.0 

2.1 

15 

50 

0 

4.5 

5.6 

5.3 

4.7 

5.3 

5.3 



0.25 

7.6 

5.7 

5.0 

6.8 

6.5 

6.8 



0.5 

5.2 

2.8 

2.1 

8.6 

6.7 

5.3 



0.75 

18.3 

1.5 

0.3 

20.0 

15.9 

4.6 


100 

0 

4.3 

4.4 

4.9 

4.2 

4.2 

4.6 



0.25 

4.7 

3.4 

3.9 

5.9 

4.3 

4.9 



0.5 

7.2 

3.8 

3.0 

9.0 

8.4 

6.5 



0.75 

40.6 

10.0 

1.8 

42.0 

36.9 

23.0 


200 

0 

4.3 

3.8 

5.4 

4.2 

4.3 

4.3 



0.25 

6.4 

5.5 

5.2 

7.7 

5.5 

6.5 



0.5 

9.3 

7.4 

6.3 

14.1 

15.7 

14.7 



0.75 

75.2 

56.5 

36.8 

79.2 

79.3 

71.5 


Table 2: Rejection rates of 1-Lq in % estimated from 1000 bivariate samples of size n generated 
under ^ 0,0 where ?^o,c and T-Li^m are defined in (1.3) and (4.1), respectively, such that, 

for each sample, the first \n9\ (resp. last n — \n9\) observations were drawn from a c.d.f. whose 
copula is the Gumbel-Hougaard, whose first margin is GEV with parameters ^ = 20, a = 10 
and 7 = 0.25) (resp. /r = 20 + d/r, a = 10 and 7 = 0.25), and whose second margin is standard 
normal. The value of the parameter of the Gumbel-Hougaard copula is set through its one-to- 
one relationship with Kendall’s tau r. 
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Figure 3: Rejection rates of the test based on ^ (o) against 0 e {0.2,0.3,..., 0.8} estimated 
from 1000 bivariate samples of size n e {100,200}, such that, for each sample, the first }ns*J 
(resp. last n — }ns*J) observations are generated from a GumbeRHougaard copula with param¬ 
eter 2 (resp. 3). The dashed line marks the corresponding estimated rejection rate of the test 
based on Sn,A- 


them to reject "Rq.c unless "Ro.m holds. Additional changes in the dispersion or scale parameter 
of the first margin might even increase the phenomenon. 

Empirical levels of the test based on ^ A consequence of Proposition 4.2 is that the 
test based on 5® ^ will hold its level asymptotically under one abrupt marginal change only, such 
as those considered in the previous experiment. To evaluate the corresponding finite-sample 
behavior, we considered again the setting of Table 1. Indeed, because of the rank-based nature 
of the test based on 5® samples generated under T-Lq can equivalently be regarded as generated 

from "Ro.c From the last two columns of Table 1, we see that the test based on 

holds its level equally well as the test based on Sn,A- The test based on S^''^ is however slightly 
too liberal for n = 50, although the agreement of its empirical levels with the 5% nominal level 
improves as n increases. 

Empirical power of the test based on As a last experiment, we investigated the 

influence of the value 6 on the power of the test based on 5"® Figure 3 displays the rejection 
rates of the test based on ^ against 6 e {0.2,0.3,..., 0.8} estimated from 1000 bivariate 
samples of size n 6 {100,200}, such that, for each sample, the first }ns*J (resp. last n — }ns*J) 
observations are generated from a Gumbel-Hougaard copula with parameter 2 (resp. 3). The 
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Pseudo-observations 



Figure 4: Left: Annual maximal peak flows and volumes of discharges measured in Streckewalde, 
Germany. Right: Corresponding pseudo-observations computed using (3.2). 


values 0.25, 0.5 and 0.75 were considered for s*. As one can see, the rejections rates are not too 
much affected by the value of 9. In addition, the power of the test based on S® ^ remains overall 
reasonably close to that of the test based on Sn,A- From a practical perspective, the latter result 
suggests that, under 77o,m in (1-2), the somehow “non optimal” use of the test based on 5®^ 
instead that based on Sn,A does not incur a large power loss, if any. As a consequence, if one 
hesitates about which of ?7o,m in (1.2) or in (4.1) holds, it seems safer to use the test 

based on 5® ^ as, should be actually true, the latter test is more likely to hold its level by 
construction, and should 'Ho,m be true, the power loss, if any, should not be too large. 


6 Illustrations 

To illustrate the proposed tests, we consider two hydrological data sets. The hrst one consists 
of n = 86 bivariate annual maxima measured between 1921 and 2011 (five years of data are 
missing) at a station located on the river Prefinitz in Streckewalde, Germany. The variables of 
interest are Q, the annual maximal peak flow (in m^/s), and V, the annual maximal volume of 
discharge (in 10® m^). Their observations are represented in Figure 4. The joint distribution of 
Q and V is of strong interest to hydrologists as it can be used to assess the risk of catastrophic 
flood levels. For a recent case study, we refer to Mitkova and Halmova (2014). 

Because we are dealing with bivariate block maxima, it is natural to assume that the data 
arise from one or more bivariate extreme-value distributions. The aim of our analysis is to 
test for possible changes in the dependence between Q and V that might have occurred during 
the long period of observation. An additional element to be taken into account here is that a 
dam was built on the river Prefinitz in 1973 (which corresponds to the 48th observation) a few 
kilometers upstream from the measurement station. We make the hypothesis that, if there are 
changes in the two components series, then, they are unique and they occurred simultaneously 
after observation 48 due to the construction of the dam. In other words, we assume that either 
'Ho,m in (1-2) holds or in (4.1) with 9 = 48/86 holds. In the former case, it is natural to use 
the test for change-point detection based on Sn,A in (3.7), while in the latter case, the extension 
based on ^ in (4.5) with 9 = 48/86 should be preferred. As mentioned in the previous 
section (see Figure 3 and the related discussion), using the test based on 5® ^ for some value 
of 9 when 'Ho.m in (1-2) actually holds does not seem to result in a strong power loss, if any. For 
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Figure 5: Left: Peak flows in rrfi/s of 164 summer flood events simultaneously measured at 
gauges in Aue and Niederschlema, Germany. Right: Corresponding pseudo-observations com¬ 
puted using (3.2). 


that reason, we carried out the test based on S'® ^ with 9 = 48/86. The resulting approximate 
p-value of 0.068, obtained from B = 10000 multiplier bootstrap replicates, indicates that there 
is some weak evidence of change in the dependence between Q and V. Interestingly enough, 
the maximum of the test statistic was not obtained for observation 48 but for observation 32 
corresponding to year 1953. 

The second data set consists of peak flows (in m?/s) simultaneously measured at two neigh¬ 
boring gauges for n = 164 physically independent summer flood events. The two gauges are 
located in Germany, in Aue and Niederschlema, respectively, and the corresponding measure¬ 
ments will thus be denoted by Qa and Qn, respectively. The observations, chronologically 
ordered, span the period 1929-2011 and are represented in Figure 5. An event was classified as 
a flood, if each peak flow exceeded the smallest annual maximal peak flow measured between 
1929 and 2011 in Aue and Niederschlema, respectively. The period of each flood event was iden¬ 
tified by hand and only the largest value (peak flow) was included in the data set. Hence, by 
construction, the observations are formed subject to a block maximal procedure, with possibly 
slightly differing block sizes for each of the flood events. It therefore seems sensible to assume 
that the data-generating distribution(s) are extreme-value distributions. 

There were two reasons why only summer events were included in the analysis. First, typical 
winter floods are produced from melting snow, whereas summer floods are due to short but heavy 
rainfalls. These very different physical mechanisms lead to different peak flow distributions. 
Second, very high peak flows, which are of particular interest, almost exclusively occur during 
the summer time. The joint distribution of peak flows is of interest, for instance, to evaluate 
the efficiency of water reservoirs (Schulte and Schumann, 2015). 

The aim of our analysis is to assess whether the dependence between Qa and Qn changed dur¬ 
ing the long observation period. As for the previous illustration, it might be important to take 
into account the fact that dams where constructed on the river Mulde and one of its tributary 
upstream of the two gauges Aue and Niederschlema. A first dam, called Schonheiderhammer, 
was put in service in 1980 (which corresponds to observation 108) and a second dam, named 
Eibenstock, was put into service in 1982. As previously, we make the hypothesis that, if there 
are changes in the two components series, then, they are unique and they occurred simultane¬ 
ously after observation 108 due to the construction of the dams. Following the same reasoning 
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as for the first illustration, we apply the test based on 5® ^ in (4.5) with 0 = 108/164 and obtain 
an approximate p-value of 0.195 based on B = 10000 multiplier bootstrap replicates. Hence, 
there is no evidence for a change in the dependence between Qa and Qn- 


A Proofs of Propositions 3.1, 3.2 and 4.1 


Propositions 3.1 and 3.2 will be corollaries of a more general result. Recall that A = {(s,s') e 
[0,1]: s ^ s'} and let An(s, s') = ([ns'J — [nsj)/n for (s, s') e A. Also, consider the process 

An{s,s',t) = v^An(s,s'){i^„^j+i,^„yj(t) - A(t)}, (s,s',t) £ A X [0,1], (A.l) 

where is defined in (3.5), and note that 

An{0,l,t) = Vn{Ann{t)-A{t)}, t e [0,1] (A.2) 

is the process of interest in Proposition 3.1. 

Theorem A.l. Under the conditions of Proposition 3.2, in the normed space (£°°(A x [0,1]), || • 
Iloo); we have A„ Ac, where 


Ac(s,s',t) =-{1 + A(t)}^ r Cc{s,s',y^ \y^)dy. 

Jo 

Proof. Since Jq l(m ^ y)dy = (1 — m) for 0 ^ m ^ 1, we can write, for 1 ^ k ^ i ^ n, 

SkAt) = 1 - f Ti ^ y} dy 

1 ^ 

^ ^ “ J £-k + i ^ ^^ ^ 

= 1 - f Ck:i{y^~\y^)dy, 

Jo 

where Ck-i denotes the empirical copula, see (3.11). Similarly, for t e [0,1], 

S(t) = 1 - {1 + A(t)}-^ = 1-C /W dy = 1 - f' C(yi-*, y*) dy. 

Jo Jo 

Therefore, introducing the notation 

Bn{s,s',t)= r Cins\+i:[ns'\{y^~\y^)dy X f ^(y^"*,y*) dy, 

JO Jo 

we obtain, for any (s, s') e A such that [ns] < [ns'J, 

'^[nsJ + l:[ns'J ‘^'(t) 1 


(A.3) 


Anis,s',t) = ^/nXnis,s') 

_ y/nXn{s ,s') 
Pn{s, S 


1 - *5[nsj + l:[ns'j(i) ^ 


f C{y^-\A)-C^nsl + l:[nsAy"-\A) 
A) Jo 

- „ / ^ f C„(s,s',y^"^y‘) dy, 

Bn{s,s',t) Jo 


dy 


(A.4) 
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where, for any (s, s', u,v) e A x [0,1]^ 


Cnis,s',U,v) = ^/n\nis,s') \^Cins\ + ^ns'v) - C{u,v)^ 

denotes the two-sided sequential empirical copula process studied in Bucher and Kojadinovic 
(2014). Since A is continuously differentiable on (0,1), we have, from Example 5.3 in Segers 
(2012), that the first-order partial derivatives Ci and C 2 of C exist and are continuous on 
(0,1) X [0,1] and [0,1] x (0,1), respectively. Hence, from Theorem 3.4 in Bucher and Kojadinovic 
(2014), 

sup |Cn(s,s',rt,u) - Cn(s,s',tt,u)| = op(l), (A.5) 

(s,s^,u,h)gA X [0,1]^ 

where, for any (s, s', u,v) e A x [0,1]^, 

C„(s, t, u, v) = ]B„(s, s', n, v) - C'i(u, u)]B„(s, s', u, 1) - C 2 {u, u)]B„(s, s', 1, u), (A.6) 

with the convention that Ci{x, •) = C' 2 (-,x) = 0 if x e {0,1}, and where 

^ [ns'l 

Bn{s,s,u,v) = -^ ^ {l{Ui^U,Vi^v) -C{u,v 

with {Ui,Vi) = {F{Xi),G{Yi)), i = l,...,n. From (A.5), the fact that ]B„(0, •,•,•) Be in 
(£”([0,1]^), II • Iloo) under ?^o,c (see, e.g., van der Vaart and Wellner, 1996, Theorem 2.12.1), the 
fact that I Cl I ^ 1 and IC 2 I ^ 1 and the continuous mapping theorem, it follows that the process 

A„(s,s',t) =-{1-t A(t)}2 X f Cn{s,s',y^~\y^)dy (A.7) 

Jo 

weakly converges in (£°°(A x [0,1]), || • ||oo) to Ae in (A.3). The theorem is thus proved if we 
show that 


sup |An(s,s',t) - A„(s,s',t)| = Op(l). 
(s,s',0eAx[0,l] 


(A.8) 


For that purpose, let us first show that, for any (s, s') e A with [nsj < [ns'J and any t e [0,1], 


Bn{s,s',t) ^ 

lb 

First of all, we have, for any u,v e [3/4,1] and any 1 ^ k < i ^ n, 

e 

£-k + l 


(A.9) 


1 ^ ^ 

Ck:e{u,v) ^ Cfc:£(3/4,3/4) = ^ l{Uk:e,i ^ ^ 3/4}. 


i=k 


The number of tjk:£,i is {£ — k + 1), and of those, exactly [3(£ — k + 2)/4j do not exceed 3/4. The 
same is true for the Hence, for at least 

[3(^ -k + 2)/4j - {(£ - fc + 1) - [3(^ -k + 2)/4j} 


of the pairs, both pseudo-observations do not exceed 3/4. As a consequence. 


Cfc,,(3/4,3/4) ^ 


2[?,{l-k + 2)/A\-{l-k + l) 
l-k + l 


3^ 

2J 


k + 2 
k 1 




1 

2‘ 
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In particular, we have Ck-.eiu^ *;?/*) ^ 1/2 for all y such that min(y^ *)?/*) ^ 3/4. Since 
min(?/^“*, y*) ^ y, we get that 



(A.IO) 


Together with the fact that Jg C{y^ ^iV^) dy = {1 + A(t)} ^ ^ 1/2, the last display implies the 
bound in (A.9). 

In order to prove (A.8), notice first that, by (A.4), (A.7), the triangle inequality and (A.9), 


sup |A„(s,s',t) 
(s,s',t)eAx[0,l] 


A„(s,s',t)| ^ 20 sup |C„(s,s',u,u)|, 

X [0,1]^ 


both sides of the inequality being zero if the suprema are restricted to (s, s') e A such that 
[nsj = [ns'J. Next, fix e, 77 > 0. Using the fact that that vanishes when s = s' and that Cn 
is asymptotically uniformly equicontinuous in probability, there exists 5 6 (0,1) such that, for 
all sufficiently large n. 


sup |A„(s,s',t) — A„(s,s',t)| > e 

(s,s^t)GAx[0,l] 


^ P "I 20 sup |C„(s, s', u, u)| > e ^ < V- 

(s,s' ,u,v)eAx[0,l]'^ 


The proof of (A.8) will thus be complete if we show that 

sup |An(s,s',t) - A„(s,s',t)| = op(l), 

(s,s',£)gAx[0,1] 

— s'^S 

which, in view of (A.9), would be an immediate consequence of the fact that 

sup \Bn{s, s',t) - {1 + A{t)}~‘^\ = Op(l). 

(s,s',£)6Ax [0,1] 
s' —s>(5 

Using again the identity {1 + C{y^~^,y^) dy, the last display follows from the fact 

that 


sup 

(s,s',£)gAx [0,1] 
s' —s5:6 


f \y^)dy-{l + A{t)} ^ 

JO 


= n-V2 


1 

sup -A 

(s,s',t)GAx[0,l] ) 

s' —s>5 


r Cn(s,s',y^ ^y‘) 

Jo 


dy 


= Op(l), 


which completes the proof. 


□ 


Proof of Proposition 3.1. The proposition immediately follows from (A.2) and Theorem A.l. 

□ 

Proof of Proposition 3.2. The assertion is a mere consequence of the fact that, under T-Lq, 
Dn(s,t) = An(s, l)An(0,S,t) - An(0, s)A„(s, l,t), (s, t) 6 [0, 1]^, (A.H) 

Theorem A.l, and the continuous mapping theorem. □ 
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For the proof of Proposition 4.1, another lemma, generalizing the identity in (A. 5), is needed. 

Lemma A.2. Under 'Hi^m ^o,c and if A is continuously differentiable on (0,1), 

sup |Cn(s,s',w,u) - C„(s, s', w,n)| = Op(l), 

(s,s',u,tj)G(An[0,6]^) X [0,1]^ 

sup \Cnis,S,U,v) -Cnis,S,U,v)\ = Op(l). 

(s,s',u,ti)s(Ari[6,l]2) X [0,1]^ 


Proof. Both suprema are measurable, and they are equal in distribution to the same suprema 
calculated under 'Ho,m T~lo,c- The assertions then follow from (A. 5). □ 

Proof of Proposition 4.1. Let 

Kis,s,t) = \/nAn(s,s'){if„^j_^i,^„^,j(t) -A{t)}, (s,s',t) 6 A x [0,1], 

where is defined in (4.3). We shall first show that 


and that 


sup IA^(s, l,t) - A„(s, 9, t) -An{9,l,t)\ = op(l) 

(s,t)s[O,0]x[O,l] 


sup |A^(0,s,t) - A„(0,6»,t) - A„(6',s,t)| = 

(s,i)s[6»,l]x[0,l] 


(A.12) 

(A.13) 

(A.14) 


where A„ is defined in (A.7). Proceeding as in (A.4), for (s,t) e [0,0] x [0,1], we have 

- S(t) 


= \/nA„{s, 1) 




_ \l,Cnis>0,y^ \y^) + Cni0.l,y^ \y^)dy 

Using the fact that 3/2 ^ {1 ^ ~ {1 + tl(t)} ^ 2 for all t e [0,1], the supremum on the 

left of (A.13) is smaller than 2/i^„ x where 


/i,„ = sup 

(s,t)e[0,e]x[0,l] 


f Cn(s,0,?/^ *,?/*)+ C„(0,l,y^ \y^) 
Jo 


dy 


and 


h,n= sup {1 - 5'[nsJ + l:n(0} -{l-'S'(i)} 

(s,0s[0,6»]x[0,l] 


-1 


From Lemma A.2, the weak convergence of Cn Cc in (f°°(A x [0,1]^), || • ||oo) under ?7o,c, 
and the continuous mapping theorem, = Op(l). Concerning l 2 ^m by the dehnition of Sf,.^ 
in (4.2), we have that 


sup 

(s,0e[0,6»]x[0,l] 


+ - S{t) = Il,n X n sup {A„(s,l)} ^ = op(l) 


sg[O,0] 


Hence, {{s,t) ^ ^ ™ (^“([Of^] x [0,1]), || • Hoo), which, from the 

continuous mapping theorem, implies that l2,n = op(l). This completes the proof of (A.13). 
The proof of (A.14) is similar. 
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To finish the proof, notice first that, under ?^o,c 

®n(s> t) = ^nis, 1)A® (0, s, t) - An(0, s)A® (s, 1, t), (s, t) 6 [0, l]^ 

where A^ is defined in (A. 12). Next, let 

]D)^(s, t) = A„(s, 1){A„(0, s A 9,t) + An{s a 0, s, t)} 

- An(0, s){A„(s, s V 6,t) + An{s V 9, (s, t) e [0,1]^, 

where a and v denote the minimum and maximum operators, respectively, 

A„(s, s',t) =-{1 + X r Cn{s,s',y^~\y^)dy, (A.15) 

Jo 


and Cn is defined in (A.6), and let 

ID)^(s, t) = A„(s, 1){A„(0, s A 9,t) + An{s a 9, s, t)} 

- A„(0, s){A„(s, s V 9,t) + An{s v 0, 1, t)}, (s, t) e [0,1]^. 

The desired result shall then follow from the fact that 


sup |D^(s,t)-ID)^(s,t)| = op(l). 

(A.16) 

(s,t)e[0,l]2 


sup |ID)^(s,t)-]D)^(s,t)| = op(l). 

(A.17) 

(s,t)G[0,l]2 


D^^Dc in (£”([0,1]2 ),||.||oo). 

(A.18) 


where Dc is given in (3.9). To show that (A. 16) holds, it suffices to restrict the supremum 
in (A.16) successively to {s,t) e [0,0] x [0,1] and {s,t) e [0,1] x [0,1] and use the triangle 
inequality, (A. 8), and (A.13) and (A. 14), respectively. The fact that (A.17) holds is obtained 
from the triangle inequality and Lemma A.2. Finally, (A. 18) follows from the fact that, for any 
0 ^ s ^ s' ^ s" ^ 1, An{s, s", •) = An{s, s', •) + A„(s', s", •), the weak convergence of C„ under 
?7o,c and the continuous mapping theorem. □ 


B Proofs of Propositions 3.3 and 4.2 

Just as for the non-bootstrap results in Propositions 3.1 and 3.2, Propositions 3.3 and 4.2 
can be conveniently proved using appropriate two-sided sequential processes. For (s,s',n, u) e 
A X [0,1]^ and b = 1,..., B, let 



u,v) 



[ns'J 

2=[nsJ + l 


and 



s', u, v) 


1 




l’^'5'J ft f . 

Ci 'I (^[nsj-I-1:[ns'J,i ^ i^[nsJ + l:[ns'J,i ^ 

i=[ns\-\-l 


(B.l) 


(B.2) 
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(B.3) 


Next, for (s, s', u,v) e A x [0,1]^ and 6 = 1,..., let 


C^^^(s, s', u, v) = lB^^^(s, s', u, v) — Ci{u, V 


s,s',y^ \l) - C2{u,v)M\^>{s,s',l,y^) 


and 


C^^(s, s', u, v) = s', u, v) — Ci(u, v)B^^(s, s', u, 1) — C 2 {u, r;)]B^^^(s, s', 1, t>). 


Furthermore, from (2.2), for {y,t) e (0,1) x [0,1], we obtain that 

Ci{y^-\y') = {^(t) -a'(t)} 

and 

C2{y^~\y') = {A{t) + (1 -t)A'(t)}y^W-*, (B.4) 

where A' is extended by continuity at 0 and 1. Indeed, from the fact that max(t, 1—t) ^ A[t) ^ 1 
for all t e [0,1], we have that A'{t) e [—1,1] for all t e (0,1), and, from the convexity of A on 
[0,1], we have that A' is increasing on (0,1). In addition, we adopt the usual convention that 
0*^ = 1. Finally, for 6 = 1,... ,B and (s, s', t) e A x [0,1], let 

=-{1+ X f C^n\s,s,y^~\y^)dy. (B.5) 

Jo 

and 

=-{1+ X f C^^\s,s',y^~\y^)dy. (B.6) 

Jo 

Lemma B.l. Under T-Lo,c and if A is continuously differentiable on (0,1), 

in (^“(A X [0,1]), |[ • lloo)^”^^; where A„ is defined in (A.15), Aq is defined in (A.3) and 
A^\ ..., a[?^ are independent copies of Ac- 


Proof. From Theorem 2.1 of Bucher and Kojadinovic (2014), the fact that \Ci\ ^ 1 and |C'2| ^ 1, 
the continuous mapping theorem and (A. 5), we have that, under ?^o,c, 

(c„,cli),...,c(f)) ^ (cc,cg\...,c(f)), 

in (^“(A X [0,1]^), II • lloo)^'*’^, where C„ and Cc are defined in (A.6) and Proposition 3.2, 
respectively, and c[)\...,C^^ are independent copies of Cc- The desired follows from the 
continuous mapping theorem. □ 

Lemma B.2. Under %q and if A is continuously differentiable on (0,1), for any b = 1,... ,B, 

sup |A®(s,s',t) - A®(s,s',t)| = op(l). 

(s,s',4)sAx[0,l] 


Proof. From the proof of Proposition 4.3 of Biicher et al. (2014) (see the term (B.3)), we have 
that, for 6 = 1,... ,B, 

sup |]B^'')(s,s',n,u)-]B®(s,s',u,'(;)| = op(l), 

(s,s^,w,hj)GAx [0,1]^ 
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where and B^^ are defined in (3.10) and (B.l). The desired result will follow if we show 
that 

sup |B^^)(s,s',n,u)-B^^)(s,s',u,u)| = op(l), (B.7) 

(s,S^,ti,tj)GAx [0,1]^ 

where B^f^ is defined in (B.2). The supremum on the left of the previous display is smaller than 


sup \C[ns\ + l-.[ns'\iu,v) - C{u,v)\ X SUp 

(s,s^u,i')gAx[0,1]^ (s,s^)gA 


L y .0^) 

n A 

2=[nsJ + l 


Using the fact that the first supremum is bounded by 2, the asymptotic uniform equicontinuity 
in probability of the process (s, s') i—>- (^y Donsker’s theorem) which vanishes 

when s = s', and the weak convergence Cc in (t'°°(A x [0,1]^), || • loo), we obtain that 

the latter display is op(l). Hence, (B.7) holds. □ 

For (s, s', t) e A X [0,1] and for 6 = 1,..., H, let 


^^nHs,s',t) = -{1 + Ai.,n{t)f X ^ y + 


(B.8) 


2 =[nsJ + l 


where ui[nsj+i:[ns'j,i is defined in (3.13). Furthermore, for 1 ^ A: ^ ^ n and (y, t) e (0,1) x [0,1], 

let 


Cl,k-J,n(l/^ ^l/^) = l^Ak:£(t) - (1 i)^ ( b . 9 ) 

C2,k:ln(y'-', y^) = {hl{t) + (1 - t)iUJA)} (B.IO) 


where A'f,.^^ is defined in (3.14) and with the convention that Ci^k:e,nix, •) = C 2 ^k-i,n{'-, x) = 0 if 
X e {0,1}. Finally, for (s, s', y,t) e A x [0,1]^ and 6 = 1,..., H, let 



{s,s',y^ *,/) 


S) S , y tV ) ^l,[ns\ + l-.[ns'\,n{y : 

^2,[ns\ + l-.[ns' \,n{y : 


s,s',y^ *,1) 
s,s', l,y‘). 


Lemma B.3. Let the pseudo-observations be either calculated as in (3.4) or in (4.4). For 
{s, s', t) e A X [0,1] and b = 1,..., B, 


'S^ A) — “{1 + Ai:n{t)} 


' r cw(s,s', 

Jo 


.1-t 


y )dy. 


(B.ll) 


Proof. From the definitions of dn, bn, Cn and given in Section 3, we have 

f Cifi(s,s',y^"^y*)dy = f B^^i(s, s', y^"*, y*) dy 

Jo Jo 

- aL«6j+i:Ln5'j(A) f (s, s', y^~\ 1) dy 

Jo 

- C[nsJ + l:[ns'j(A) f (s, s', 1, y^) dy 

Jo 

= (J'l - h- h){s,s',t), 
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where /i, I 2 and I 3 are defined in an obvious manner. Observe that we can write 


[ns'J 


^ ^ {^[nsJ + l:[ns'J,i(^) ^ 2/} 

2^[nsJ + l 


1 


[ns'J 


[ns'J — [nsj 


2 1 {m 

[nsJ+l:[ns'J,jr (0 ^ v} 


j=[ns\+l 


Then, since Jq l(m ^ y)dy = {1 — m) for 0 ^ m ^ 1, we obtain that 

^ [ns'J 

) y) j — ^ j {^[nsJ + l:[ns'J (^) ^[nsj + l:[ns'j,i(0} • 


j=[nsj + l 


Similarly, since Jg y^ ^l(m ^ y) dy = (1 — m^)lb, 


h{s,s,y) = 


- .. [ns'J 

®[nsJ + l:[ns'J T1 1 ^(b) t— 


and 


h{s,s',y) = 


^[nsj + l:[ns'j(^) ^ j=[nsj + l 

*^[nsJ + l:[ns'J ItJ 1 ^(b) f- 


Ci {^[nsJ + l:[ns'J (0 f^[nsJ + l:[ns'J,i(^)} 


{^[nsJ+l:[ns'J (^) ^[nsJ + l:[ns'J,i(^)} • 


'^[nsj + l:[ns'j(^) i=[nsj + l 

Hence, (B.ll) is proved. □ 

Proof of Proposition 3.3. For (s, t) e [0,1]^ and b = 1,..., H, let 

= An(s,l)A^^)(0,s,t) - An(0,s)A®(s,l,t), 
where A® is defined in (B.5), and let 

= An(s,l)A^'’)(0,S,t) - An(0,s)A^'’)(s,l,t), 
where An^ is defined in (B.6). From Lemma B.2, we then immediately obtain that 

sup |lD)^'’^(s,t) -]D)®(s,t)| = op(l). 

(s,t)G[0,l]2 

The latter combined with Lemma B.l, the continuous mapping theorem, (A. 5), (A. 8) and (A.11) 
gives 

(d„,dW,...,d(,^)) ^ (Dc,og\...,o[f^) 

in (l'“([0,1]^), II • where Dc is defined in (3.9) and are independent copies 

of Dc. From the definitions in (3.12) and (B.8), we further have that, for {s,t) e [0,1]^, 

®nHs,t) = An(s,l)A^'’)(0,S,t) - An(0,s)A®(s,l,t). 


Hence, to complete the proof, it remains to show sup(g_^)g|-o ^2 |ID)n^(s,t) — Dn''('S,t)| = op(l), 
which is implied by 




sup |A^'’)(s,sT) - A^^)(s,s',t)| = op(l). 
(s,s',t)GAx[0,l] 


(B.12) 
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Having in mind the fact that Ai-n converges uniformly in probability to H as a consequence of 
Proposition 3.1, (B.7) and the fact ICil ^ 1 and IC2I ^ 1, to prove (B.12), it suffices to show 
that 


sup 

(s,s^)eA 


{f^l,[nsJ + l:[ns'J,n(y lU ) , 


s,s',y^ *,1) = op(l) 


(B.13) 


and 


sup 

(s,s')eA 

(!/,t)6[0,l]2 


{C2Xns\ + l-.\ns'\,n{y \y)-C2{y \ , s\ I, y^) 


op(l)) 


(B.14) 


where Ci,[nsj+i:[ns'j,n(y^ C 2 Xns\+i\ns'ln{y^ *,!/*), Ci{y^ \y^) and C 2 {y^ ^y*) are given 

in (B.9), (B.IO), (B.3) and (B.4), respectively. 

From (A.10), it can be verihed that, for any n ^ 1 and any t e [0,1], |Ai:„(t)| ^ 7. Since, 
by definition, in (3.14) is also uniformly bounded, {y,t) >—>■ Ci^k:e,niy^~^,y^) and {y,t) 

C 2 ,k-i,n{y^~^,y^) are uniformly bounded. 

To prove (B.13), we can then proceed as in the proof of Proposition 4.3 of Bucher et al. (2014) 
(see the terms (B.4) and (B.5)). Using the asymptotic uniform equicontinuity in probability of 
which vanishes when s = s', and the fact that Ci and its estimator are uniformly bounded, 
it remains to show that, for any <5 e (0,1), 


sup 

(s,s^,y,t)€Ax [0,1]^ 
s' — s'^8 


Cl\ns\ + l-.\ns'\,ni.y iV ) ^l(2/ ) 2/ ) 


Op(l). 


(B.15) 


The previous result is implied by the fact that 


sup 

{s,s',i)6Ax[0,l] 
s' — s'^5 


2^[nsJ + l:[ns'J ^(^) 


sup 

(s,s',t)eAx[0.1] 

s' — s'^5 


An{s, S , t) 

^/nXnis,s') 


op(l)) 


where is defined in (A.l), and an analogue result for „ dehned in (3.14). The 

latter can be seen as follows: for ^[„sj+i [ns'j n ™ (3-15), we have 


sup l^[nsJ+l:[n.'J,nW -^'Wl ^ SUp 

s,s',t)eAx[hn,l-hn] tG[h„,l — hn] 


A{t + hn) - A{t - hn) 


2hr. 


- A'it) 


+ sup 

(s,s',t)EAx [hji,l — hji\ 


S -\- hn) ^ •) ^ hn) 


2hn\fn\n{s, s') 


. (B.16) 


Since A' , extended by continuity, is (uniformly) continuous on [0,1] (see the discussion be¬ 
low (B.4)), by the mean value theorem, the first term on the right converges to zero. The 
second term on the right is smaller than 


sup |A„(s, s', t) — A„(s, s', t')| X 

{s,s',t,t')eAx[0,l]2 

t-t'|sS2hn, 


1 

sup “——7 
(s,s0eA ‘^hn^nyS^ S ) 

s' —s'^5 


op(l), 


by asymptotic uniform equicontinuity in probability of A„. 

Hence, (B.15) holds and, thus, so does (B.13). The proof of (B.14) is similar. This completes 
the proof of (B.12) and, therefore, of the proposition. □ 
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Proof of Proposition 4.2. For b = 1, ..., B, let be the analogue of A® in (B.8) based 
on the adapted pseudo-observations defined in (4.4). Then, by definition of we have 

= A„(s,l)A^’(')(0,s,t)-A„(0,s)A^’(')(s,l,t), (s, t) 6 [0, l]^. 

Next, let 

= An(s, 1){A^'’H0>'S A 0,t) +A^n\s A 9,S,t)} 

- A„(0,s){A®(s,s V 6i,t) -h A^^^(s v 6», l,t)}, {s,t) e [0, l]^ 

where An^ is defined in (B.6), and let us first show that 

sup |D^’(')(s,t)-ID)^'(')(s,t)| =op(l). (B.17) 

(s,t)s[0,l]2 

To prove (B.17), we shall show that 


sup 

(s,t)G[O,0]x[O,l] 


and 


sup 

(s,i)e[0,l]x[O,l] 


+A„(O,s){AW(s,0,t)+AW(0,l,t)}| =op(l), (B.18) 

- An(s, l){A^^)(0,6',t) -h A®((9,S,t)} 

+A„(0,s)AW(s,l,t) =op(l). (B.19) 


We start with the proof of (B.18). Under (B.12) continues to hold if the supremum 

is restricted to (s, s', t) e (An [0, 6*]^) x [0,1]. This can be seen by essentially the same arguments 
as in the proof of Lemma A.2. Therefore, (B.18) will hold if 


sup 

(s,i)e[0,e]x[0,l] 


A^’(') (s, 1, t) - A® (s, 9, t) - aW {9, l,t) = op(l), 


that is, if 


sup 

S6[O,0] 

(u,«) 6 [ 0 , 1]2 


B^’(')(s,l,u,u)-BW(s,0,u,u)-BW(0,l,n,u) = op(l), (B.20) 


and, having in addition (B.3), (B.4), (B.9) and (B.IO) in mind, if 

sup \^[ns\ + l-.nit) - Mt)\ = Op(l), 

(s,t)s[o,e]x[o,i] 

sup \^[ns\ + l:n,nit) “ = Op(l), 

(s,Le[0,6»]x[0,l] 


(B.21) 

(B.22) 


where is the analogue of in (3.14) defined from the adapted pseudo-observations 

in (4.4). The supremum on the left of (B.20) is smaller than 


sup 

(s,u,7;)g[O,0] X [0,1]^ 


[n9\ - [nsj ^ ^ ^ 

n-[ns\ ‘^M+i:Kl(uw) 


n — {n9\ 
n — [resj 


C[ne\ + l-.n{u,v) - C{U,V) 


X sup 
se[O,0] 


1 ” 

7 ^ , 2 , 


(b) 


^ sup 

(s,u,t))s[O,0] X [0,1]2 


i=[nsj + l 

\Cn{s,9,U,v) + Cn{9, 1,U,U)| 

y/nXn{s, 1 ) 


X Op(l) = op(l) 
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by Lemma A.2 and the weak convergence of Cn under T-Lq^c The supremum on the left of (B.21) 
is smaller than 

sup {A„(s, 

(s,t)s[0,6»]x[0,l] 

where A® is defined in (A.12), and is op(l) because of (A.13), (A. 7), Lemma A.2 and the weak 
convergence of under 7^o,c- The proof of (B.22) is based on a decomposition similar to that 
used in (B.16) and relies again on (A.13). Hence, (B.18) holds. The proof of (B.19) is similar. 

Using arguments of the same nature as those employed in the proof of Lemma A.2, we obtain 
the following extension of Lemma B.2 under 

sup |A^'’^(s,s',t) - A^^)(s,s',t)| =op(l), 

(s,s',t)E{(An[0,6»]2)u(An[e,l]2)}x[0,l] 

By the triangular inequality, this implies that 

sup (s, t) I = op(l), (B.23) 

(s,t)s[0,l]2 

where 

= A„(s, 1){A^'’)(0,S A e,t) A 9,S,t)} 

- An(0,s){A^^^(s,s V 0,t) + V 9,l,t)}, {s,t) e [0,1]^. 

and An'^ is defined in (B.5). The desired result is then a consequence of Lemma B.l, the 
continuous mapping theorem, (A.16), (A.17), (B.17) and (B.23). □ 


C Test statistic and multiplier bootstrap for d ^ 2 


In Sections 3 and 4, we restricted ourselves to the case d = 2. Results for arbitrary dimension d ^ 
2 can be established at the cost of a more complex notation but without significant additional 
mathematical difficulties. We give the main steps of the generalization hereafter. Let X = 
(Ai,..., Xfi) be a random vector with c.d.f. and extreme-value copula of the form (2.1) and (2.2), 
respectively, and suppose that A is continuously differentiable on the interior of Sd-i with partial 
derivatives Aj{t) = dA{t)/dtj, j = 2,... ,d. With the notation Uj = Fj{Xj), j = 1,... ,d, and 
ti = ti{t) = I — ^ ^ •Sd-i, we have, just as for d = 2, 

A{t) = 5(t)/{l — S{t)} and S{t) = E ( max 

with the convention that = g for all u e (0,1). 

Let Xi, i = 1,... ,n, he independent copies of X and let Uk:£,i = ■ ■ ■ ■, Uk-i,id) be d- 

variate generalizations of the “subsample” pseudo-observations in (3.4). We define a CUSUM- 
type process on [0,1] x by 

, Inshn—Inshfi ,, ^ .fo 

I®n('S,t) = -- |^l:[nsj(*) “ ^[nsj + l:n(^) j ) 

where, for 1 ^ A; ^ ^ n, Akx{t) = Sk-£{t)/{l — Sk-iit)}, and 



Sk-A't) 


1 

--;- > max 

i — k + 1 l-gjjgd 

i=k 


fjAh 

^k:e,ij 
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with the convention that j = 0 if /c > £. 

Let us introduce some additional notation. For any y e [0,1] and t e Sd-i, we define y* 
to be the vector e [0,1]'^ with the convention that 0^ = 1. Furthermore, for any 

u e [ 0 , l]'^ and any j = 1,... ,d, denotes the vector of [ 0 , l]'^ whose components are all 
equal to one, except the jth which is equal to Uj. 

Proposition C.l. Suppose that all of the above conditions are met. Then, in the normed 
space .^“([0,1] X Dc, where 

= {1 +X r sCc{s,l,y*) - {I-s)Cc{^,s,y*)dy. 

Jo 

Here, Cc denotes a centered Gaussian process on A x [0,1]*^ defined through 

d 

Cc(s,s',u) = {Bee'S', u) -Be(.s,u)} - ^ Cj(u){Bc(s',u‘^^^) -Bc(s,u'^^^)}, 

i=i 

where Be is a tight, centered Gaussian process on (f“([0, || • |[oo) with covariance kernel 

given by 

E{Be(.s,it)Be(s', u')} = (s a s'){C'(mi a u'l,... ,Ud c\ u'd) — C{u)C{u')} , 
and Cj, j = 1,..., d, denotes the jth first-order partial derivative of C. 


The proof is almost identical to that of Proposition 3.2. For a corresponding bootstrap 
approximation of the limit De, let i = 1 ,... ,n, b = 1 ,... ,5, be i.i.d. standard normal 
multipliers. Furthermore, from (2.2), we have that, for any y e (0,1) and t e 5 ^- 1 , 


Cjiy^) = 


y 






yA{t)-tAA{t)+Afit)-^‘\,^^tyA^ 


j = 1 , 

, j = 2,..., d. 


The above quantities can be estimated consistently by plugging in subsample estimators of A 
and Aj, j = 1 ,..., d, respectively, namely A^-i and 

~ ^2h~ , J = 2,..., d, 

with tBhnBj = {t 2 ,... ,tj-i,tj±hn,tj+i,... ,td) and a sequence hn i 0 such that inf^^i hn^/n > 
0 (boundary effects can be dealt with by generalizing the approach adopted below (3.15)). Then, 
analogously to the bivariate case, we define 

[ns] u\ 

/ J '*^[nsj + l:n,i(*) ^3/2 ^ ^l:[nsjd(i) j", 


(s, t) — (1 + Ai-riit)}'^ X 
where, for 1 ^ k ^ ^ n, 


n 


3/2 


2=[nsJ + l 


2 = 1 


Wk-.i,i{t) = mk-.e{t) - rhk,(.,i{t) + ^ 

1=1 


ifdk\i^ij ( 2 ) Uk:£,j ( 2 ) 


bk:£,j ( 2 ) 

with fhkj and Uk:£j denoting the arithmetic mean over i = k,... ,£ of 
rhk:e,i{t) = max and Uk:£,ij{t) = 
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and where 


+ ■^j,k-l,n{i') ^j/=2^j' ^j' :k:£,n{t) 1 j ~ ■ ■ • i d-i 

^kl,j{t^ ~ “1“ 1 tj. 

Test statistics and corresponding multiplier bootstrap replicates can be defined analogously to 
Section 3, as functionals of and b = 1,..., B, respectively. In addition, generalizations 
adapted to known breaks in the margins can be obtained by computing pseudo-observations 
from the subsamples determined by the marginal change-points, as explained in Section 4. We 
omit the details for the sake of brevity. 
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